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ABSTRACT 

The Cartwheel galaxy harbours more Ultra-Luminous X-ray sources (ULXs) 
than any other galaxy observed so far, and as such it is a particularly inter- 
esting target to study them. In this paper we analyse the three Chandra 
observations of the brightest ULX (NIO) in the Cartwheel galaxy, in light 
of current theoretical models suggested to explain such still elusive objects. 
For each model we derive the relevant spectral parameters. Based on self- 
consistency arguments we can interpret NIO as an accreting binary system 
powered by a ^ 100 Mq black hole. A young supernova strongly interacting 
with its surroundings is a likely alternative, that can be discarded only with 
the evidence of a flux increase from future observations. 
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1 INTRODUCTION 

The Ultra-Luminous X-ray sources (ULXs) are a class of very bright {Lx — 10^^—10^^ erg s~^), 
point-like X-ray sources detected off the nuclei in several galaxies (Fabbiano 2006). Al- 
though they were discovered more than 20 years ago (by the Einstein X-ray satellite: 
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Long & van Speybroeck 1983; Fabbiano 1989), their nature remains unclear. If the ULXs 
are powered by accretion, their luminosity exceeds the Eddington limit for a stellar-sized 
object, sometimes by a factor ~ 1000. This has led some authors (Colbert & Mushotzky 
1999) to postulate that the ULXs are black holes with masses 10^ - 10"^ Mq, intermediate 
between the black holes expected from the final evolution of massive stars and the super- 
massive black holes of 10^ — 10^ Mq powering the Active Galactic Nuclei. Possible formation 
scenarios of these intermediate mass black holes (IMBH) are discussed in Miller & Hamilton 
(2002) (formation in globular clusters), and Portegies Zwart ct al. (2004) (in young super- 
massive star clusters). An alternative view places the ULXs in the more familiar realm 
of stellar-sized black holes. According to King et al. (2001), the ULXs are black holes of 
M, ~ 10 Mq, accreting mass from a disc, at a rate above their Eddington limit. In this 
regime, the inner accretion disc is thick, and may coUimate the outgoing radiation. The 
anisotropy leads the observer to overestimate the luminosity (and the mass of the hole) by 
a large factor, up to 10 — 100. Alternatively, some models suggest that the ULXs may ac- 
tually radiate above the Eddington limit. They are slim disc models (Ebisawa ct al. 2003), 
photon bubble dominated discs (Begelman 2002, 2006) and two-phase radiatively efficient 
discs (Socrates & Davis 2006). A combination of beaming and super-Eddington accretion is 
also possible (Poutanen et al. 2007; King 2008, 2009). 

More recently, a third possibility has been emerging, namely that the ULXs are black 
holes of 30 — 90 Mq produced by the final evolution of massive stars in a metal poor environ- 
ment (Zampieri & Roberts 2009). In such an environment the black hole's progenitor does 
not loose much of its original mass by the action of line-driven winds, and might collapse 
directly into a black hole without exploding as a supernova, producing more massive black 
holes than in a metal rich environment. 

The accretion-powered models do not exhaust the possible scenarios able to explain 
the nature of the ULXs. In a sample of 154 ULXs observed with Chandra, Swartz ct al. 
(2004) argue that about 20% of them may be well described by a thermal model, consistent 
with young supernova remnants strongly interacting with the surrounding environment. In 
this interpretation, the Eddington limit is clearly not an issue, and indeed the ULXs have 
X-ray luminosities comparable to those of the brightest supernovae {Lx ~ 10"^^ erg s~\ 
Immler & Lcwin 2003). 

The Cartwheel galaxy is a peculiar nursery of ULXs. This galaxy underwent a collision 
with a smaller galaxy about 10® yr ago (Higdon 1996; Mapelli et al. 2008). This episode not 
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only conferred the galaxy its peculiar shape, but also triggered a massive star formation in its 
ring. The Cartwheel harbours the largest number of ULXs than any other galaxy: 15 sources 
are more luminous than 10^^ erg s~^ (Wolter & Trinchieri 2004), and at least some of them 
are known to be variable (Crivellari et al. 2009). In this paper we address the properties of 
the ULX labelled NIO by Wolter & Trinchieri (2004). In the Chandra observation of 2001 
it was the brightest ULX in the Cartwheel (and also among the brightest known ULXs), 
but a couple of subsequent XMM-Newton observations taken in 2004 and 2005 showed 
that its luminosity is fading (Wolter et al. 2006). Two additional Chandra observations 
were performed in 2008 to follow the variability pattern of this ULX. We present here the 
Chandra observations taken in 2008 and compare them with that of 2001. The ACIS angular 
resolution is instrumental to reduce to a minimum the possible contamination of NIO from 
the surrounding diffuse gas and the neighbour sources crowding the Cartwheel ring. In our 
analysis we do not use the XMM-Newton observations, since the analysis of the variability 
with XMM-Newton data was already published by Wolter ct al. (2006). In addition, the 
use of XMM-Newton data would require specific modelling of the contamination of the 
neighbouring sources and the Cartwheel's ring due to the relatively large PSF of XMM- 
Newton, incrementing the uncertainties in the derived spectral parameters. 

Our main goal is to compare the spectral parameters of NIO to the theoretical models 
put forward to explain the ULXs. In particular, we shall discuss several accretion models and 
the supernova model. Although NIO is a bright source, its large distance {D = 122 Mpc, 
Wolter & Trinchieri 2004 and references therein) limits the number of collected photons, 
preventing us from rejecting or confirming a spectral model on purely statistical grounds. 

For this reason, we shall assess the likelihood of any model more in the light of its own 
self-consistency than of its statistical evidence. 

The outline of the paper is the following. In section 2 we describe the Chandra data and 
their preparation. In Section 3 we present a simple spectral model of NIO, which will be a 
thread for a more detailed analysis. The accretion and supernova models are presented in 
section 4 and 5, respectively. Finally, we discuss and summarise our results in section 6. 

2 DATA PREPARATION 

Chandra observed the Cartwheel once in 2001 and twice in 2008. All the observations were 
carried out with the back-illiminated (BI) chip ACIS-S3. This detector was operated in 
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FAINT mode in the first observation (no. 2019) and in VFAINT mode in the remaining 
observations (nos. 9531 and 9807). 

New level=2 event files have been re-created with CIAO-4.1.2 in order to work with a 
homogeneous data set. The script wavdetect was used to detect the sources within the CCD 
field of view. No observation was seriously contaminated by flares; to check this we have 
first extracted a light curve from the source-free areas in the spectral band 0.3 — 10 keV. 
The Chandra script Icsigma-clip used this curve to excise the periods with a count rate 
more than 3 sigmas above the average. This correction excludes but a small fraction of the 
observations, as shown in the summary Table 1. 

The spectra of NIO were extracted from a circle of radius R = 2.7", and the background 
spectrum from a source-free region near the source. The chosen extraction circle includes 
virtually all the photons from the source. In all our analysis we fit the spectra in the window 
0.3 — 7 keV, where the signal to noise ratio is higher. In this window the first spectrum 
(observation no. 2019) contains 400 photons, the second one (observation no. 9531) 173 
photons and the third (observation no. 9807) only 59 photons before background subtraction. 
There are two contributions to the background: instrumental and due to the diffuse gas of 
the galaxy ring and to unresolved low luminosity sources. Both of them give a negligible 
contamination on account of the small extraction area: we estimate that the instrumental 
background contributes for less than 2 photon in each observation, and 3 — 5 photons come 
from the gas of Cartwheel galaxy's ring. 

On account of the relatively poor statistics of our data, we would not like to loose spectral 
resolution by binning the data to the standard minimum value of ~ 25 counts/bin required 
by a consistent use of the statistics (see e.g. Cash 1979). For this reason we choose to 
bin the spectra with a minimum of 5 counts/bin and analyse them with the implementation 
of the Cash statistics (called cstat) provided in the version 12.5.1 of XSPEC. The Cash 
statistics assumes that the counts are distributed according to a Poisson law, which does not 
allow to subtract the background from the spectra (Cash 1979). In our case the background 
contribution is small, and we simply neglect it. 

The cstat statistics differs from the canonical Cash statistics because the value of cstat 
provides a goodness of fit (similar to the statistics), if each spectral bin contains at least 

5 counts. The fit is considered acceptable if the value of the "reduced" cstat (i.e. the ratio 
of ctstat to the number of degrees of freedom of the model) is close to one (Arnaud 1996 
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^). We also check that the distribution of residuals in any spectral model is not skewed and 
does not show peculiar features. 

All the statistical uncertainties of the best-fitting parameters are quoted at 90% confid- 
ence level. The X-ray luminosities of the source are consistently calculated in the spectral 
band 0.5 — 10 keV, and they are corrected for absorption. 

3 THE POWER LAW MODEL 

The count rate between the three observations varies. In order to avoid complications in 
comparing count rates of observations taken in FAINT and VFAINT modes, we evaluate 
this variability with a simple absorbed power law model (wabs * powerlaw). Despite its 
simplicity, it is informative of some basic parameters useful to our further analysis. Separate 
fits of each observation return comparable values of the absorption column uh and the 
spectral photon index F. For this reason we tie uh and F of each data set to a common 
value and fit them simultaneously. Table 2 summarises the results of the fit, and Figure 1 
shows the model and its residuals. The value of the C-statistics is C/dof = 100.8/101. The 
best fitting absorption column is uh = 'i-7^io x 10^^ cm~^, and the spectral photon index 
is F = 1.88^0 24, both consistent with those presented by Woltcr & Trinchieri (2004). The 
value of Uh is higher than Galactic, but consistent both with other X-ray measures in the 
Cartwheel and with the intrinsic absorption in the optical band (see Wolter & Trinchieri 
2004). 

The intrinsic luminosities in the spectral band 0.5 — 10 keV is 1.2 X 10^^ erg s-\ 8.6 x 
Xq40 gj^g g-1 g^j^j^ 2.8 X 10^° erg s~^ for the first, second and third observation respectively. 
Compared to the Eddington luminosity of a 1 Mq collapsed object, these luminosities seem 
to imply a black hole of > 200 Mq. It is perhaps worth to remark that the Eddington 
luminosity is bolometric, while in the present paper we refer to the luminosity in the X-ray 
band. This is an approximation of the bolometric luminosity within a factor > 2. 

Formally, a power law model provides an adequate description of the Chandra data sets, 
and more refined models are not statistically required. Nevertheless, we would like to know 
whether these observations may somehow constrain some of the theoretical models suggested 
for the ULXs. The following sections aim at this goal. 



See also the XSPEC manual page http://heasarc.gsfc.nasa.gov/docs/xanadu/xspoc/manual/XSappendixCash.html. 
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4 THE ACCRETION MODELS 

In the framework of an accretion scenario, the engine of NIO is a black hole accreting mass 
from a donor star through a disc. In this section we present the results of the spectral analysis 
of three accretion models: i) a multi-colour disc (diskbb in the language of XSPEC), ii) a 
model of a disc around a maximally rotating Kerr black hole (kerrd in XSPEC) and 
finally iii) a Kawaguchi (2003) slim disc model ^. All models are acceptable on statistical 
grounds and they all return similar values of the C-stat. For an easier comparison, we have 
summarised the results of these models in Table 3. We also shall discuss the hyperaccretion 
model suggested by King et al. (2001), but since no spectral model is available for this, our 
discussion will be mainly qualitative. 

4.1 The multicolour disc 

The first model we consider is a multicolour disc (MCD) modified by the interstellar absorp- 
tion. The two free parameters of the model are the effective temperature of the inner rim 
Tin and the apparent inner disc radius R^^, related to the actual inner radius -Rdisc through 
the relation (Kubota et al. 1998) 



where ^ ^ 0.412 is a correction factor that takes into account that the maximum temperat- 
ure Tin does not occur exactly at the inner edge of the disc (Kubota et al. 1998). For typical 
parameters, the disc optical opacity is dominated by electron scattering. The electron scat- 
tering may also Comptonise the emerging spectrum, and this effect is taken into account by 
introducing the so called hardening factor /, i.e. the ratio / = Tcoi/Tcs between the colour 
temperature Tcoi and the effective temperature Teg. The value of / is insensitive to the disc 
viscosity and the mass of the accretor, but is sensitive to the mass accretion rate. For a mass 
accretion rate close to the Eddington limit the value of / lies in the range / ^ 1.7 — 2.0 
(Shimura & Takahara 1995). Kawaguchi (2003) argues that if the accretion rate largely ex- 
ceeds the Eddington limit [Mc^/L-Edd > 100), the hardening factor is somewhat larger, in 
the range / ~ 2.3 — 6.5. Introducing the radiation efficiency 



-Rdisc — ^ f Ri 



(1) 



r] = L/Mc^, 



(2) 



^ Slim disc models arc not available in the main release of the package, but only as tabular models to be downloaded separately 
from the URL http://heasarc.gsfc.nasa.gov/docs/xanadu/xspec/models/slimdlsk.htnil. 
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this limit translates into 

L>100r7LEdd. (3) 
For a disc accreting on a Schwarzschild black hole r} ~ 0.06 (Frank et al. 2002), so Kawaguchi's 
hardening parameter applies if > 6LEdd- 

The normalisation of the model only involves geometrical parameters, so we link the 
normalisation of the three spectra to a common value, and leave the disc temperatures free 
to vary. Since the absorption columns inferred from each spectrum are consistent, we also link 
them to a common value. The interstellar absorption column isun = l.S^Qg x 10^^ cm~^. The 
disc inner temperature is = 1.33lo^7 f*^^ ^^'^ ^^^^ observation; Tin = l-2ll;o}g keV for 
the second observation and Tin = 0-89toii keV for the third observation. The normalisation 
of the model is i^' = 7.4I33 x 10~^, and is related to the apparent inner radius i?in through 
the expression 

= 1.22 X 10^ km ( ' ( ^ (cos^)"'/', (4) 

where 6 is the inclination of the disc with respect to the line of sight = is face-on), and 
D is the distance to the source. Inserting the appropriate value of K we find 

i?disc ^ 3.3tli X 10^ e f/icosOy/^ km. (5) 
If the inner edge of the disc is located at the last stable orbit of a Schwarzschild black hole, 
this implies 

M. ~ 61+}^ ( —^—rj^. (6) 

° V2.O J \0A1 J (cose)i/2 ^ ^ 

The resulting disc luminosities are Li = 4.1 x 10^° erg s~^, L2 = 2.7 x 10^° erg s~^ and L3 = 
8.0 X 10^^ erg s~^, for the first, second and third observation, respectively. All these values 
exceed the Eddington luminosity TEdd — 7.9 x 10'^^ erg s~^ for a black hole of M, ~ 60 Mq 
by a factor L/L-^dd ^ 5, consistent with our choice of the hardening parameter / = 2. 

A long-standing problem of the application of MCD models to the ULXs is that they 
tend to overestimate the disc's inner temperature Tin. Indeed, the luminosity of a MCD is 
(Kubota et al. 1998) 

^boi = 4vr asB (^)' (7) 

(Makishima et al. 2000). If the temperature is too high (for the same luminosity and the 
other parameters), the disc radius -Rdisc is underestimated, so is the black hole's mass. 
A possible explanation for this effect is the presence of a hot corona embedding the disc 
(Stobbart et al. 2006): a single MCD model applied to such system would return spuriously 
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high disc temperatures. We investigate how the presence of a hot corona would affect our 
estimate (6) of the black hole mass. We describe NIO with a Comptonised disc emission, 
where the temperature of the Compton seed photons is set equal to the disc inner tem- 
perature. On account of our rather low statistics, the temperature of the hot corona is not 
well constrained, and it has been fixed to ksTc = 50 keV, while its optical depth has been 
left free to vary. The best-fitting value of the Compton thickness is r ^ 0.5lo:2: this results 
strongly depends on the temperature Tc, with cooler coronae returning larger optical depths. 
The temperature of the disc, however, is not sensitive to the precise values of Tq and r. As 
noted by Stobbart ct al. (2006), a Compton-thick (r > 1) corona obscures the inner region 
of the disc, and in this case the MCD parameters cannot provide reliable estimates of the 
black hole's mass (see also Gladstone et al. 2009). The disc temperature of the three obser- 
vations are respectively Tin = O.SS^JJ;}^ keV, Tin = 0.36lJ^:J^ keV and Tin = 0.09l[5:|^| keV. 
The normalisation of the disc component is K = 7.6^l^g^ x 10"'^ (in XSPEC units), yielding 



Using the above disc inner temperatures we infer the luminosity of the disc component to be 
Li = 9.0 X 10^^ erg s~^ (first observation), T2 = 1.5 x 10^^ erg s~^ (second observation) and 
L3 = 5 X 10'^^ erg s~^ (last observation). The ratio between Li and the Eddington luminosity 
of a black hole of M, ~ 141 Mq is ~ 0.5, which justifies our choice / = 1.7 of the hardening 
parameter. As expected (see Stobbart ct al. 2006), a hot corona lowers our estimate of disc's 
temperature and increases the black hole's mass with respect to a simple MCD model. 

4.2 The Kerr disc 

The second model we consider is the spectrum emitted by a disc orbiting a maximally 
rotating Kerr black hole. We fix the distance of the Cartwheel to its known value, and the 
hardening parameter to / = 1.7: the quality of our fit is independent of this choice. The 
inclination of the disc with respect to the line of sight has been fixed to ^ = 0°; the value 
of 6 affects the determination of M,, as we shall discuss at the end of this section. We also 
set the radius of the accretion disc to the last stable orbit of a maximally rotating Kerr BH 
(^ 1.235 G M,/c^), and the disc outer radius to its default value (the model is insensitive to 
this parameter). The only free parameters of the model (apart from the hydrogen column 
density) are the mass M, of the central object and its mass accretion rate M. In the fitting 




(8) 
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procedure we link all the parameters for the three observations except the mass accretion 
rates. 

The fit returns the absorption column nfj = 2.2lQg x 10^^ cm~^, and the mass of the 
central object results 

M. = 92.81^^;^ Mq. (9) 

The luminosity of the source is Li = 4.3 x 10"^*^ in the first observation, L2 = 2.9 x 10'^° in 
the second, and L3 = 8.5 x 10'^^ in the last. The Eddington luminosity associated to the 
mass (9) is LEdd = 1-2 x 10^'' erg s~^. These values must be checked for their consistency with 
the adopted hardening parameter. A Kerr maximally rotating black hole has an accretion 
efficiency rj ~ 0.4 (Frank ct al. 2002, see Equation 2). I<Cawaguchi's larger values of / then 
apply if L > 40LEdd- In all our observations, the Kerr disc model returns L/L^dd ^ 3.5, so 
the adopted hardening parameter is consistent with the inferred M, and the relatively low 
mass accretion regime. 

Before ending this section we need to address the effect of the disc inclination 6 on the 
estimate of the mass M,. The value of 6 does not affect the quality of the fit, but the mass 
is sensitive to it, as shown in Figure 2 (see also Hui & Krolik 2008). Low inclination angles 
imply lower M,, which (for / = 1.7) attains a minimum value of ~ 77 Mq for 6 = 20°. If 
the disc is seen almost edge-on, on the other hand, M, grows up to M, ~ 10"^ Mq. This 
trend of M, with the viewing angle is also found in the multicolour discs considered in the 
previous section (in particular, see Equation 6). 

To summarise, the Kerr disc model is consistent with a black hole of M, c:^ 90 Mq, 
similar (within the uncertainties) to what inferred from the Comptonised multicolour disc 
model. Higher masses are not excluded, if the disc is highly tilted with respect to the line 
of sight. Also a higher hardening factor would increase M,, but higher values of / are not 
required by the inferred accretion regime. 

4.3 Slim Discs 

The multicolour disc model assumes that the disc radiates locally as a black body. If the ac- 
cretion rate largely exceeds the Eddington limit, this hypothesis is not appropriate. In such 
an accretion regime the heat content of the disc is trapped and dragged to the black hole's 
horizon before it is radiated. For this reason, slim discs are radiatively inefficient, and their 
properties significantly differ from those of the standard thin discs (see e.g. Abramowicz ct al. 
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1988; Frank et al. 2002). In this section we investigate whether the data support the possib- 
ility that NIO is powered by a slim disc accreting above its Eddington limit. 

In XSPEC two shm disc models are available: we adopted the model based on the 
calculation by Kawaguchi (2003), with the mass of the central object, the accretion rate and 
the viscosity a as free parameters. The spectrum is calculated taking into account the local 
Comptonisation and the relativistic effects. Kawaguchi 's model has been applied to fit the 
XMM-Newton spectra of some ULXs (see e.g. Okajima et al. 2006; Vierdayanti et al. 2006), 
returning black hole masses of few tens of Mq accreting at super-Eddington regimes. In 
our fit we have linked the interstellar absorption column, the hole mass and the a parameter 
between the observations. The mass accretion rate is independent for each observation. The 
fitting procedure is unable to constrain the black hole's mass. The confidence regions are 
open for M, > 80 Mq (Figure 3), and the model is consistent with any value of M, above 
this limit. The absorption column is uh = 3.6^Qg x 10^^ cm~^, while the viscosity parameter 
and the BH mass are a = 0.56_o.37 and M, = 495_34o Mq respectively. The mass accretion 
rates for the three observations are Mi = 39.4_23, M2 = 21.41^2^0, M3 = S.ltJ^g' units 
of MEdd = L-Edd/c"^- Assuming the conversion factor MEdd = 1-3 x 10^^ g s~^(M,/100 Mq) 
used by Kawaguchi (2003), the best fitting accretion rates are Mi = 2.5 x 10^^ g s~\ 
M2 = 1.4 X 10^^ g s"\ M3 = 3.3 X 10^° g s"^ It is not possible to set upper limits to the 
the parameters a, M, and Mi, since the error calculation hits the limits of their tabulated 
values (a = 1, M. = 1000 Mq and M = 1000 MEdd, respectively). 

One caveat is in order about the consistency of the application of slim disc models to 
NIO. The radial temperature profile of slim discs varies as T oc R^^/'^ instead of T oc 
^-3/4 characteristic of thin multicolour discs (Watarai et al. 2000). The radial dependence 
of Tin in NIO may be checked with a multicolour disc (diskpbb in XSPEC) where the 
temperature index p (defined by T{R) oc R~^) is a free parameter. The index p ~ 0.5 has 
been found in some ULX spectra successfully fitted with a slim disc model (Okajima et al. 
2006; Vierdayanti et al. 2006). In our case the model wabs*diskpbb is comparable to the 
others and since its best-fitting parameters are very similar to those derived from a standard 
multicolour disc, we do not present them. Our best fitting value is p ~ 0.74l;Q'^g close to the 
standard value 0.75, but slightly inconsistent with the slim disc value p = 0.5. 
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4.4 The "Hyperaccretion" model 

The last accretion model we consider for NIO is the so-called "hyperaccretion" scenario 
(King et al. 2001; King 2008, 2009). This model explains the high luminosity of the ULXs 
with a combination of a high accretion rate (close to the Eddington limit) on a stellar mass 
black hole, and a mechanical beaming due to the accretion stream itself. 

The apparent X-ray luminosity (Equation (10) in King 2009) of the source results 



L ~ 2.2 X 10^^ erg s"^ 



M,\ M 



2 r 



MqJ \MEdd 



In 



M 



Ml 



Edd 



(10) 



where M is the mass transfer rate from the donor star and Medd = L-^aa/v^^^ is the 
Eddington-limited mass transfer, dependent on the efficiency r] of the X-ray production 
in the accretion process, defined by Equation (2). 

One of the most relevant features of this model is that the beaming factor scales as M~^. 
This causes the bolometric luminosity to scale with the disc's temperature as 

L oc Tr\ (11) 

opposite to the standard correlation L oc T^. 

The prediction (11) has been found consistent with the L — Tin relation of the soft excess 
observed in some ULXs (Feng & Kaaret 2007; Kajava & Poutanen 2009; Gladstone et al. 
2009). We assume that the count rate of NIO correlates with its intrinsic luminosity in all 
our observations. In principle, it is certainly true that a strong spectral change may alter 
(even reverse) this correlation, but such a variation of the spectrum of NIO is not supported 
by our data. Under this working hypothesis, we measure a positive correlation between the 
inner temperature and the luminosity, at odds with the prediction of the hyperaccretion 
model. The hyperaccretion scenario probably is not the correct explanation for the nature 
of NIO. 



5 THE SUPERNOVA MODEL 

As discussed in the introduction, accretion power is not the only viable explanation for 
the ULXs' engine. Alternative models are possible, and young SNe interacting with their 
surroundings may explain an appreciable fraction of the known ULXs (Swartz ct al. 2004). In 
this section we investigate the SN model for NIO, aiming at deriving the spectral parameters. 

The X-ray emission from supernovae may start up to about one year after the explosion, 
and it is powered by the interaction of the ejecta with the circumstellar medium (GSM, see 
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e.g. Immler & Lewin 2003). Shock fronts form as the ejecta expand through the CSM. The 
leading (or "forward") shock heats the CSM up to 10^ — 10^° K, and gradually an inner (or 
"reverse") shock wave starts to propagate. The density behind the reverse shock is 5 — 10 
times higher than that behind the forward shock, and the temperature is lower, around 
10^ — 10*^ K (Chevalier & Fransson 1994). For this reason, the soft X-ray emission (below 
< 5 keV) is dominated by the reverse shock, and is well modelled by thermal models, while 
the hard emission (say, above 10 keV) is dominated by the forward shock. 

The luminosities of X-ray supernovae lie in the range 10'^'' — 10^^ erg s~^, with the so- 
called SN Iln supernovae at the bright end of the interval. These are core-collapse super- 
novae, owing the "n" in their name to the presence of several narrow emission lines in their 
optical spectrum (Turatto 2003). These lines are most likely due to the interaction of the 
ejecta with a dense, slow wind emitted by the SN progenitor. Other supernovae explode 
in thinner environments, resulting in dimmer X-ray luminosties. For this reason, if NIO is 
actually a young SN, it is most likely a SN Iln. 

We fit the spectra of NIO with the thermal models apec and nei, corrected for inter- 
stellar absorption. The apec model assumes that the emitting plasma is in full collisional 
equilibrium: there is a dynamic balance between the collisional ionisation and the electronic 
recombination. The nei model does not assume collisional equilibrium, allowing a mismatch 
between the collisional ionisation and the (longer) recombination process. This mismatch 
may prevail in shocked plasmas, where the electrons were energised but did not yet settle to 
collisional equilibrium with the ions (see e.g. Dopita & Sutherland 2003). The parameters of 
the nei model are very similar to those derived for the apec model (the fit is quite insensitive 
to the ionisation age), so we present only the apec results. 

In our model (wabs*apec) the absorption column and the metal abundances of the 
three observations are linked together. Since there is no compelling statistical evidence of 
a variation of the temperature, we also link together the temperatures of the three ob- 
servations. The results of this model are presented in Table 4. The statistics of the fit is 
C/dof = 95.9/100, and the values of the best-fitting parameters are riu = 2.^^^^ x 10^^ cm^, 
ksT = keV, Ni = 4.3^5^2 ^ ^0^^ (normalisation of the first observation), N2 = 

3.1^09 X 10~^ (normalisation of the second observation), and = 1.0lg2 ^ 10~^ (nor- 
malisation of the third observation). The units of these quantities are (for sources at low 
redshift) 10~^^ / dVUf, uh/Ati D"^, where D is the distance to the source, rig and uh are the 
electron and proton densities of the emitting plasma, and the integral (known as emission 
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integral) is extended over the volume occupied by the source. The emission integrals are 
consistent with the size (few 10^^ cm) and the density (~ 10"^^ — 10^^^ g cm~^) expected 
from the reverse shock of a young supernova remnant few years after the explosion (see e.g. 
Chevalier & Fransson 2003). 

The unabsorbed light curve of NIO, calculated according to the power law model, is 
plotted in Figure 4 (see also Woltcr ct al. 2006). The evolution of the observed flux cannot 
be fitted by a power law, as observed in other X-ray supernovae (1986J, Houck ct al. 1998; 
1988Z, Schlegel & Petre 2006; 1998S, Pooley et al. 2002). This is not an argument against 
the SN nature of NIO, though, since a simple power-law decline of the X-ray luminosity 
has not been reported in other confirmed SN Iln. SN 1995N, for instance, shows a complex 
light curve, possibly resulting from the interaction of the ejecta with a complex circumstellar 
medium (Zampieri ct al. 2005). The sudden decrease of the X-ray luminosity of NIO in 2008 
might show that the ejecta have passed the dense region of the pre-supernova wind, and 
have finally reached a low density outer environment. 

SN Iln are usually found in high metal environments (1986J, Houck ct al. 1998; 1998S, 
Pooley et al. 2002). Indeed, high metal abundances favour the blowing of the line-driven 
strong pre-supernova winds required to set up the dense circumstellar environment required 
to power a stronger X-ray emission. The metallicity of the Cartwheel is low as measured 
from the HII regions (Fosbury & Hawardcn 1977); the X-ray measure is poorly constrained, 
and consistent with zero. Therefore the SN progenitor had to be quite massive to blow dense 
enough winds without the aid of a high metallicity. 



6 DISCUSSION AND SUMMARY 

In this paper we considered the spectral modelling of the ULX NIO, located in the Cartwheel 
galaxy in order to assess the nature of this source. Although the source is intrinsically very 
bright, the spectra have few counts on account of the large distance (122 Mpc). For this 
reason, we are unable to reject (or confirm) a model on purely statistical grounds. We first 
discussed several accretion models (multicolour disc around a Schwarzschild or Kerr black 
hole, slim disc, hyperaccretion). All models indicate a black hole of ~ 100 M0, at the high 
end of the mass distribution of a black hole generated by stellar evolution, possibly in a metal 
poor environment. This result is consistent with the conclusions of Mapelli et al. (2009) 
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and Zampieri & Roberts (2009) that stellar evolution in a metal depleted environment may 
produce black holes of 30 — 90 Mq. 

The theoretical studies of the evolution of very massive stars have been rekindled by 
the discovery of the ULXs, with the result that black holes of ~ 100 Mq are more com- 
mon than earlier works suggested. Isolated solar-abundance stars may have masses up to 
150 Mq, but in dense environments they may coalesce to form stars as massive as 1000 Mq 
(Yungelson ct al. 2008). Recent studies on the mass loss from hot, massive stars (Puis ct al. 
2008, and references therein) have shown that the occurrence of weak and/or clumpy winds 
may reduce the mass loss up by a factor 10, thus allowing the star to retain a higher fraction 
of its original mass at the end of its life. Also lower metallicities entail weaker wind mass 
losses (Heger et al. 2003). Zero metallicity (Population III) stars may be significantly more 
massive than solar-abundance stars (Ohkubo ct al. 2006). All these factors determine the 
mass of the final black hole, that may range from ~ 70 Mq for a solar-abundance progen- 
itor (Yungelson et al. 2008), up to 500 Mq for low metallicity stars (Ohkubo ct al. 2006). 
Observations indicate that the low metallicity scenario is more likely for NIO: the Cartwheel 
galaxy is metal-poor, and although the available measurements (Fosbury & Hawarden 1977) 
are not at the exact location of NIO, they should nevertheless be representative of the aver- 
age abundance in the ring. We therefore conclude that the interpretation of NIO as a black 
hole of ~ 100 Mq suggested by the X-ray data is consistent with the theoretical predictions 
of the evolution of massive stars in a dense environment. 

In summary, if ordinary stellar evolution is the correct scenario to form the black hole, 
NIO could be an extreme High Mass X-ray Binary (HMXB). The mass accretion rate M 
on the black hole may be inferred from the luminosity L via Equation (2), and is of the 
order of M ~ 10"^ Mq yr^^. This value is comparable with the loss rate of massive (i.e., 
few tens solar masses, Frank ct al. 2002) donor stars on a thermal time scale. The accretion 
flow from the donor star to the BH is most likely to occur through a Roche lobe overflow; 
the BH capture of a wind blown by the companion is less favoured, since the small capture 
radius would require an unlikely strong mass loss from the companion. 

Can the observed decay of the luminosity of the source be explained in the framework of 
the accretion model? The longest characteristic time of an accretion disc is the so-called "vis- 
cous time" tvisc, i-e. the characteristic time taken by the disc to adapt to new conditions. For a 
standard a-disc orbiting a black hole of 100 Mq with an accretion rate M ~ 10^^ Mq yr^^, 
tvisc — 10^ — 10'^ s (Frank et al. 2002), much shorter than the time scale of the variabil- 



Chandra observations of the ULX NIO in the Cartwheel galaxy 15 

ity of NIO (few years). Therefore, the observed variabihty cannot be due to a sort of disc 
instabihty, which would affect the disc and the luminosity of the source over a time scale 
tvisc- The simplest alternative is that the variability is due to a change of the mass transfer 
rate from the donor star. This may be due to an intrinsic decay of the mass loss from the 
donor, but also to other effects occurring at the inner Lagrangian point, since the instant- 
aneous mass flow here is very sensitive to the relative sizes of the donor star and the Roche 
lobe (Frank et al. 2002). New observations of this interesting source would help to settle the 
question. 

We also explored the possibility that NIO is a young supernova strongly interacting with 
the circumstellar medium. Both the interpretation of this model and the best-fitting spectral 
parameters are consistent with this view. A new flux brightening of the source in the future 
would rule out this hypothesis. 
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Figure 1. Data and residuals for the best-fit absorbed power law for the three Chandra, observations. Top (black) obs. no. 2019, 
middle (red) no. 9531, bottom (green) no. 9807. 
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Figure 2. The best-fitting mass Af, of the model wabs * kerrd as a function of tiie inclination between the disc and the line 
of sight. The error bars represent the 90% uncertainties calculated with the free parameters. 
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Figure 3. The dependence of the fit C-statistics on the black hole mass in the slim disc model by Kawaguclii (2003). Since 
the C-statistics curve is flat for M, > 80 A/q , the best-fitting mass of the hole can be anywhere above this value. 
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Figure 4. The light curve of the source NIO. Fluxes have been calculated with a power-law model corrected for the interstellar 
absorption. The symbols refer to ROSAT data (triangle) XMM-Newton (circles) and Chandra (squares). ROSAT and XMM- 
Ncwton data are taken from Woltcr et al. (2006). 
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Table 1. Chandra, observation log of the Cartwheel galaxy. 



Observation 
identification 


Exposure time 
(before) after cleaning 
(ks) 


Count rate 
(0.3 - 7 keV) 
(10-3 cts/s) 


Mode 


Start date 

(yyyy/mm/dd) 


2019 


(76.5) 76.1 


5.3 


FAINT 


2001/05/26 


9531 


(52.8) 48.0 


3.6 


VFAINT 


2008/01/21 


9807 


(49.8) 49.5 


1.2 


VFAINT 


2008/09/09 



Table 2. Summary of the absorbed power— law model (wabs * powerlaw). 



Obs. id. 


Normalisation 
10-^ XSPEC units 


-Fb.5-10 kcV 
(10-" erg s-i cm-2) 


-^^0.5-10 kcV 
(10*" erg s-i) 


2019 
9531 
9807 


^■■^^-0.27 


6.60 
4.81 

1.56 


11.8 
8.6 
2.8 


UH = 3.67tj5-|^ X 10^1 cm-2 


r = i.88ij;-i^ 


C/dof = 100.8/101 





iimmary of the accretion disc models; MCD=wabs*diskbb, MCD+C=wabs* (diskbb+comptt) , Kerr=wabs*kerrd Slim=wabs*Kawaguchi. The first group of rows shows the best-fitting 
3r three observations separately. The second groups presents the values of the parameters common to all the observations, and the other parameters discussed in the text. Fluxes and 
ire calculated in the spectral band 0.5 — 10 keV. Kp and Kq stand for the disc and comptt normalisation, respectively. /, r) and are the hardening parameter, the accretion efficiency 
ration of the disc. The luminosites are total, so the luminosity of the MCD+C model also includes the contribution of the Compton corona. 





MCD 






MCD+C 






Kerr 






SUm 




( keV) 


Flux 
10-14 
erg s-l cm— 


Lx 
10-10 
erg s^i 


( keV) 


Flux 
10-14 
erg s-l cm— 


Lx 
lO-io 

erg s-l 


M 

10^0 g s-l 


Flux 
10-14 
erg s-l cm— ^ 


10« 

erg s^i 


M 


Flux 
10-14 
erg s-l cm-^ 


1040 
erg s-l 


1 OO+0.22 

1.00_Q j7 

0.89±?,;it 


4.63 
3.09 
0.88 


4.1 
2.7 
0.8 


5^+" " 

U.OO_Q j2 

S6+°1'5 

m+°'°^ 

u.uy_o Q4 


5.32 
3.63 
1.13 


8.6 
6.3 
2.0 


7 1+0.8 
' •-'--0.7 

4.8«:« 

1 c:+0.4 


4.84 
3.28 
0.94 


4.3 
2.9 
0.85 


39.4-23 
21.4+^2^ 
5 1+79 


11.5 
8.0 
2.3 


5.8 
4.0 
1.1 




7.4^3.3 X 10-* 


2.0 
0.06 






1 7+0.8 
-■-•'-0.7 

7.6liOg'5 X 10-3 
2.27+0% X 10-^ 
50 keV (fixed) 
5+°-4 


1.7 
0.06 






1 (fixed) 


92.8t2?:| 
1.7 
0.4 






6.72 X 10-9 (fixed) 


0.56-0.37 

495 - 340 
1.7 
0.06 






92.8/101 






88.5/99 






93.2/101 






96.7/100 
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Table 4. Summary of the absorbed APEC model wabs* apec. 



Obs. identification no. 


Normalisation 
10"^ XSPEC units 


-Fb.5-10 kcV 
(IQ-i'* erg s^i cm"2) 


i'O.5-10 kcV 
(10*0 erg s-i) 


2019 
9531 
9807 


4 3+0.9 
o 1+0.7 

■'•^-o.g 


5.55 
4.04 
1.31 


9.9 
7.2 
2.3 


UH = 2.8lo;g ^ 10^^ '^^'^ 


kgT = 5.llf;^ keV 


Z = 0+0-2 


C/dof = 95.9/100 
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